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The parsec-scale radio properties of blazars detected by the Large Area Telescope (LAT) on board 
the Fermi Gamma-ray Space Telescope have been investigated using observations with the Very Long 
Baseline Array (VLB A). Comparisons between LAT and non-LAT detected samples were made using 
contemporaneous data. In total, 232 sources were used in the LAT-detected sample. This very large, 
radio flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that 
produces strong 7-ray emission. It has been found that LAT-detected BL Lac objects are very 
similar to the non-LAT BL Lac objects in most properties, although LAT BL Lac objects may have 
longer jets. The LAT flat spectrum radio quasars (FSRQs) are significantly different from non-LAT 
FSRQs and are likely extreme members of the FSRQ population. Contemporaneous observations 
showed a strong correlation, whereas no correlation is found using archival radio data. Most of 
the differences between the LAT and non-LAT populations are related to the cores of the sources, 
indicating that the 7-ray emission may originate near the base of the jets (i.e., within a few pc of the 
central engine). There is some indication that LAT-detected sources may have larger jet opening 
angles than the non-LAT sources. Strong core polarization is significantly more common among 
the LAT sources, suggesting that 7-ray emission is related to strong, uniform magnetic fields at the 
base of the jets of the blazars. Observations of sources in two epochs indicate that core fractional 
polarization was higher when the objects were detected by the LAT. Included in our sample are 
several non-blazar AGN such as 3C84, M82, and NGC 6251. 



I. INTRODUCTION 

The Large Area Telescope (LAT; |7j) on board the 
Fermi Gamma-ray Space Telescope is a wide-field tele- 
scope covering the energy range from about 20 MeV 
to more than 300 GeV. It has been scanning the entire 
7-ray sky approximately once every three hours since 
July of 2008. The 7-ray bright sources from the Fermi 
LAT First Source Catalog (1FGL; 0) that have been 
identified with radio sources are most often associated 
with blazars (685 of 1451 sources). These blazars typ- 
ically are strong, compact radio sources which exhibit 
flat radio spectra, rapid variability, compact cores 
with one-sided parsec-scale jets, and superluminal mo- 
tion in the jets [TtJ • The dominance of blazars contin- 
ues in the Fermi LAT Second Source Catalog (2FGL; 
21]), compiled using the first two years of LAT data. 

Several very long baseline interferometry (VLBI) 
programs, such as the Monitoring Of Jets in AGN 
with VLBA Experiments (MOJAVE; [15]) observ- 
ing at 15 GHz, the Boston University program ob- 
serving at 22 and 43 GHz and the Track- 
ing AGN with Austral Milliarcsecond Interferometry 
(TANAMI; [Hf) observing at 8.4 and 22 GHz, along 



with the LAT collaboration themselves @ are in place 
to monitor the radio jets from the brightest blazars 
such as 3C273, BL Lac, etc. Our sample has a flux 
limit roughly an order of magnitude below the MO- 
JAVE survey and so allows us to probe the extensions 
of the radio core/7-ray properties down to a fainter 
population. 



II. SAMPLE DEFINITION 

Our contemporaneous radio sample is made up of 
Very Long Baseline Array (VLBA) observations made 
in 2009 and 2010, all at a frequency of 5 GHz. We 
made new observations of 90 sources from the VLBA 
Imaging and Polarimetry Survey (VIPS; [10]) and 
new 5 GHz observations of 7 sources in the MO- 
JAVE sample. We refer to this VIPS follow-up sample 
as "VIPS+". We also observed 135 non-VIPS, non- 
MOJAVE sources selected from the Fermi 1FGL as 
sources which were associated with a source in the 
Combined Radio All- Sky Targeted Eight GHz Survey 
(CRATES; @) with high (>80%) probability and had 
a flux density of at least 30 mJy in CRATES. We refer 
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to this sample of 135 new objects as "VIPS++" . The 
VIPS++ sample expands the number of LAT sources 
in our sample and includes objects outside of the orig- 
inal VIPS target RA and DEC, which were limited to 
the region covered b y th e Sloan Digital Sky Survey 
Data Release 5 (see [10] more information). Of our 
232 sources, 95 are BL Lac objects, 107 are FSRQs, 
and 30 are other types of AGN (radio galaxies, AGN 
of unknown type, and 1 starburst galaxy). Any ob- 
ject that is not a BL Lac object or FSRQ we clas- 
sify as 'AGN/Other'. The optical classifications are 
taken from the LAT First Catalog of AGN (1LAC; 
Q). See [lij for more information on this contempo- 
raneous sample. 

To build a non-LAT detected sample for compar- 
ison, we excluded all LAT sources from VIPS, leav- 
ing 1018 objects in our non-LAT sample. Of these 
1018 objects, 24 are BL Lacs, 479 are FSRQs, and 
515 are AGN/Other types (radio galaxies or AGN of 
uncertain type). The optical types for VIPS sources 
were adopted from the Candidate Gamma-Ray Blazar 
Survey (CGRaBS; [§]). It should be noted that the 
non-LAT sample is not contemporaneous with LAT 
observations. 



III. GAMMA-RAY AND RADIO FLUX 

The LAT 7-ray fluxes are reported in several energy 
bands in both 1FGL and 2FGL. To create total 7-ray 
fluxes, we combined the fluxes from 3 bands: 100-300 
MeV, 300 MeV - 1 GeV, and 1-100 GeV. The fluxes 
were added and uncertainties were added in quadra- 
ture. However, some sources had only upper limits to 
their fluxes in certain bands. If a source's reported 
fluxes in one or two bands were upper limits, we use 
1/2 the reported flux as the uncertainty in that band 
as the upper limits are given as 2-sigma results Q. 

A. GAMMA-RAY - RADIO CORRELATION 

In Figure [1] we plot the 1FGL LAT flux versus the 
total VLB A flux density at 5 GHz. Again, the LAT 
fluxes are broadband fluxes from 100 MeV to 100 GeV. 

We used the nonparametric Spearman test [23[ to 
look for correlation between the LAT flux and the to- 
tal VLBA flux density. The Spearman p s values we 
found were 0.467 for the BL Lacs and 0.510 for the 
FSRQs. The probabilities for getting the same p s val- 
ues from random sampling were 2xl0 -6 for the BL 
Lac objects and 2xl0~ 8 for the FSRQs. Therefore, 
the LAT flux correlates very strongly with the total 
VLBA flux density. The BL Lac objects do not corre- 
late quite as strongly, especially considering that we 
found a marginal correlation between total VLBA flux 
density and redshift for these objects. Still, this pro- 
vides solid evidence that objects with higher radio flux 



10' 



10 



10 







A ■ 

» i 

■ * 


• 


. . '? ' 






■ BL Lacs 

* FSRQs 

• Other - 


■ 



10 10 10 10 

Total VLBA flux density at 5 GHz (mjy) 



FIG. 1: LAT 7-ray flux (100 MeV - 100 GeV) vs. total 
VLBA radio flux density at 5 GHz. The 7-ray fluxes are 
from 1FGL and are in units of 10 -9 photons cm -2 s _1 . 
The black squares are BL Lacs, the blue triangles are FS- 
RQs, and the red circles are radio galaxies and unclassified 
objects. 



density also produce more 7-ray flux. 

In we found that there was no strong correla- 
tion between radio flux density and 7-ray flux. How- 
ever, our flux densities for those observations were 
prior to the launch of Fermi. Furthermore, it is likely 
that the radio flux density increases during episodes 
of 7-ray flaring [l^ . Our newer, larger, and more 
contemporaneous set of observations are more appro- 
priate for making these comparisons. 



B. LAT VERSUS NON-LAT RADIO FLUX 
DENSITIES 

The median total VLBA flux density at 5 GHz was 
177 mjy for the LAT BL Lacs and 221 mjy for the 
non-LAT BL Lacs. The LAT FSRQ appeared to have 
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FIG. 2: Distributions of LAT (top) and non-LAT (bottom) 
FSRQ total VLBA flux density at 5 GHz. 
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significantly higher flux densities than the non-LAT 
FSRQs, with a median of 467 mJy for LAT FSRQs 
compared to a median of only 191 mJy for non-LAT 
FSRQs. The LAT AGN/Other objects also had higher 
5 GHz flux densities than their non-LAT sources coun- 
terparts. The Kolmogorov-Smirnov (K-S) test [23j] 
probability that the LAT and non-LAT BL Lac ob- 
jects belong to the same parent population is about 
38%, although we should note that we have a very 
small sample of non-LAT BL Lac objects. The K-S 
test probability that the LAT and non-LAT FSRQs 
are related is only 3xl0 -12 . See Figure [5] for a plot 
of the FSRQ total VLBA 5 GHz flux density distri- 
butions. 
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C. GAMMA-RAY LOUDNESS 

For a measure of "7-ray loudness," we use the ratio 
of the 7-rayluminosity to the radio luminosity as de- 
scribed in [16j]. However, unlike [l6j], we do not have 
multiple observations of our sources from which to de- 
termine a median radio luminosity. Instead, we only 
use a single VLBA observation to calculate a luminos- 
ity. Due to the variable nature of these blazars, this 
may not be the best representation of the sources' ac- 
tual average luminosity. For the 7-ray luminosity, we 
used the average fluxes reported in the 1FGL and the 
2FGL. Whenever possible, we used the redshifts re- 
ported in 1LAC. For sources with no redshift in 1LAC, 
we searched the NASA/IPAC Extragalactic Database 
(NED). 

All of our sources were significantly 7-ray loud, with 
7-ray /radio luminosity ratios of at least 400 u sing 
both 1FGL and 2FGL data. It was reported in [16|] 
that the MOJAVE group saw a significant difference 
between the 7-ray loudness distributions for BL Lac 
and FSRQs. We did not find any indication of this 
in our data. The K-S test returned a 8% probability 
that the 1FGL BL Lac objects and FSRQs were drawn 
from the same parent population. For the 2FGL data, 
the K-S test result was 18%. 



IV. PARSEC-SCALE RADIO PROPERTIES 
A. CORE BRIGHTNESS TEMPERATURE 

We looked at the distributions of the core brightness 
temperatures for our LAT and non-LAT objects. The 
LAT FSRQs were significantly different from their 
non-LAT counterparts. The K-S test gave a 6 x 10 -8 
chance that the two distributions were drawn from the 
same parent population. The median core brightness 
temperatures were 6.4 x 10 10 K for the LAT FSRQs 
and 2.5 x 10 10 K for the non-LAT FSRQs, indicating 
that the LAT FSRQs are more extreme sources. The 



FIG. 3: Distributions of LAT (top) and non-LAT (bottom) 
core brightness temperatures. The BL Lac object distri- 
butions are in black on the left. The FSRQ distributions 
are in blue on the right. 



BL Lac objects did not show a significant difference 
between the LAT and non-LAT populations. See Fig- 
ure [3] for the distributions of core brightness tempera- 
tures for the BL Lac objects and FSRQs. A difference 
between the core brightness temperature distributions 
was also reported in [l2| , although they did not sepa- 
rate objects by optical type. 

We also looked for correlation between core bright- 
ness temperature and the frequency at which the syn- 
chrotron radiation is at a maximum (Vp eak ) for our 
LAT-detected objects. Whenever possible, we used 
published values of v^ eak (from Q], 0], [H, and HI). 



For those sources without published v% 



peak ' 



used the 



estimation method described in |4j, based on the av- 
erage radio-optical and optical- x-ray spectral indices 
published in [6j . If we could not find a published value 
or use the estimation method, we obtained data from 
NED and fit a parabola in log space. We only found 
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FIG. 4: Core brightness temperature versus peak syn- 
chrotron frequency for LAT-detected sources. 
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FIG. 5: Percent of sources with detected core fractional 
polarization in various observing campaigns. VIPS ob- 
servations were made during and prior to 2006. VIPS+ 
observations were made in 2009-2010, VIPS++ observa- 
tions were made in 2010. Gray indicates the sample with 
no LAT-detected sources, black indicates LAT samples. 



a significant correlation for the BL Lac objects. The 
Spearman test resulted in a ps of —0.4 and a p of 
10 . See Figure [4] for a plot of core brightness tem- 
perature versus peak synchrotron frequency. It looks 
very similar to the results presented in [16]. 



B. CORE FRACTIONAL POLARIZATION 

In order for a source to be considered "polarized" we 
required that it have a polarized flux of at least 0.3% 
of the Stokes / peak value (to avoid leakage contam- 
ination) and have at least a 5c detection (compared 
to the noise in polarized flux density). See [14j for 
more discussion on the polarization poperties of our 
sources. We found that the LAT sources are polarized 
much more often than the non-LAT sources (see Fig- 
ure EJ. The abundance of polarized sources among 
the 7-ray emitting population indicates that strong, 
uniform magnetic fields at the base of the jets are 
somehow related to the production of 7-rays. Also, 48 
of the 75 polarized sources which were in both VIPS 
and VIPS+ showed an increase in core fractional po- 
larization during LAT detection. 



C. JET OPENING ANGLE 

We measured the a mean apparent opening half an- 
gle for each source with core-jet morphology following 
the procedure described in [2J]. This method relies 
on averaging the apparent opening angles found for 
each jet component in a source. That is, the apparent 
opening angle is measured for several (more than 2) 
jet components. The final mean opening angle for the 
source is the average of all of the opening angles for 
the individual jet components. 
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FIG. 6: 2FGL 7-ray to radio luminosity ratio versus ap- 
parent jet opening angle. The 1FGL data did not show a 
correlation. 



The MOJAVE group reported in [lfj that they 
found a non-linear correlation between the apparent 
opening angle and the 7-ray loudness. We also found 
a correlation, however only when we used the 2FGL 
data. When we performed the Spearman test on 
all 34 sources with both a measured opening angle 
and a 2FGL 7-ray flux, we found a p$ of 0.57 and 
a p of 0.00044. When we break the sources up by 
type, we only see a tentative correlation for the FS- 
RQs (ps=0.58, p=0.0088), however the sample sizes 
are very small. We only had 36 objects with open- 
ing angle and 7-ray loudness (G r ) measurements in 
1FGL, and 34 in 2FGL. The MOJAVE group's sam- 
ple in [16| is well over 100 sources. See Figure H] for 
a plot of the 2FGL 7-ray loudness versus opening an- 
gle. We should also note that the MOJAVE group in 
[161 ] made their measurements of the apparent opening 
angle using mean sizes of jet components over several 
epochs, not the mean of multiple apparent opening 
angles from several components in a single observa- 
tion as we did. The fact that MOJAVE is a moni- 
toring program makes it better suited to investigating 
the apparent (and intrinsic) opening angles than our 
sample. 



D. JET CHARACTERISTICS 

We measured the change in jet position angle ("jet 
bending" , see [TTj] ) and the length of the jet (the maxi- 
mum separation between the core component and the 
jet components). We applied the K-S test to each 
of these measurements to see if the LAT and non- 
LAT distributions were different. We did not find any 
significant differences for either of these properties. 
However, we do note that the LAT BL Lac objects 
have a "long-jet" morphology (a jet with an angular 
size of at least 6 milliarseconds) more often than the 
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FIG. 7: Distributions of jet brightness temperatures for 
LAT (top) and non-LAT (bottom) FSRQs. 

non-LAT BL Lac objects (see Our measure- 

ments of the jet brightness temperature (the bright- 
ness temperature of the brightest jet component), on 
the other hand, did produce an interesting result. We 
found that the LAT FSRQs had significantly higher 
jet brightness temperatures than the non-LAT FS- 
RQs. The K-S test returned a 10~ 5 chance that the 
LAT and non-LAT distributions were drawn form the 
same parent population. The FSRQ jet brightness 
temperature distributions are shown in Figure [7] 

V. CONCLUSIONS 

Our LAT BL Lac object sample is nearly 4 times 
the size of our non-LAT BL Lac object sample. We 
only find 3 minor differences between the two: (1) 
the LAT BL Lac objects are polarized more often, 
(2) the LAT BL Lacs are more often of a "long-jet" 
morphology, and (3) the LAT BL Lac objects may 
have larger apparent jet opening angles. This leads 
us to suspect that all BL Lac objects are emitting 
7-rays. 

Our LAT and non-LAT FSRQ samples are often 



very different. The LAT FSRQs appear to be ex- 
treme sources. They tend to have high radio flux den- 
sity, high core and jet brightness temperatures, and 
(maybe) large apparent jet opening angles. The LAT 
FSRQs are also polarized much more often than the 
non-LAT FSRQs. 

The correlations found between the radio flux den- 
sity and 7-ray flux indicate that the synchrotron and 
inverse Compton radiation are related. Therefore, the 
7-rays should originate in the jets, where we know the 
synchrotron-emitting electrons are located. However, 
because most of the differences between the LAT and 
non-LAT populations are associated with the cores of 
our blazars, we suspect the 7-rays are emitted at the 
base of the jets (i.e., within a few parsecs of the cen- 
tral engine) . The prevalence of polarized cores among 
the LAT-detected sample indicates that strong, uni- 
form magnetic fields play a role in the production of 
the 7-rays. 
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